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In this paper we study the c-function of the sine-Gordon model taking explicitly into account the 
periodicity of the interaction potential. The integration of the c-function along trajectories of the 
non-perturbative renormalization group flow gives access to the central charges of the model in the 
fixed points. The results at vanishing frequency /3 2 , where the periodicity does not play a role, are 
retrieved and the independence on the cutoff regulator for small frequencies is discussed. Our findings 
show that the central charge obtained integrating the trajectories starting from the repulsive low- 
frequencies fixed points (/3 2 < 8-7 t) to the infrared limit is in good quantitative agreement with the 
expected Ac = 1 result. The behavior of the c-function in the other parts of the flow diagram is also 
discussed. Finally, we point out that also including higher harmonics in the renormalization group 
treatment at the level of local potential approximation is not sufficient to give reasonable results, even 
if the periodicity is taken into account. Rather, incorporating the wave-function renormalization (i. 
e. going beyond local potential approximation) is crucial to get sensible results even when a single 
frequency is used. 
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I. INTRODUCTION 

Statistical field theory has known an outpouring de¬ 
velopment in the last decades, with the systematic im¬ 
provement of powerful theoretical tools for the study of 
critical phenomena and phase transitions [lj. A key role 
among these methods is played by the renormalization 
group (RG) approach: RG allowed to treat statistical 
physics models studying the behavior of the transforma¬ 
tions bringing the microscopic variables into macroscopic 
ones W- Using RG one can have not only a qualitative 
picture of the phase diagram and fixed points, but also 
accurate quantitative estimates of critical properties as 
critical exponents and universal quantities - even though 
very few cases are known where RG procedure can be 
carried out exactly and the method itself offers few pos¬ 
sibility to obtain exact results. 

On the other hand, the scale invariance exhibited by 
systems at criticality may give rise to invariance under 
the larger group of conformal transformations @] locally 
acting as scale transformations [5]. The conformal group 
in d spatial dimensions (for d ^ 2) has a number of in¬ 
dependent generators equal to \{d+ l)(d + 2), while for 
d = 2 the conformal group is infinitely dimensional [5]. 
The occurrence and consequences of conformal invariance 
for 2-dimensional field theories have been deeply investi¬ 
gated and exploited to obtain a variety of exact results 
m and a systematic understanding of phase transitions 
in two dimensions. 

A bridge between conformal field theory (CFT) tech¬ 
niques and the RG description of field theories is pro¬ 
vided in two dimensions by the c-theorem. Far from fixed 
points, Zamolodchikov’s c-theorem [B can be used to get 


information on the scale-dependence of the model. In 
particular the theorem states that it is always possible 
to construct a function of the couplings, the so-called c- 
function, which monotonically decreases when evaluated 
along the trajectory of the RG flow. Furthermore, at the 
fixed points this function assumes the same value as the 
central charge of the corresponding CFT [I]. 

Although the c-theorem is by now a classical result, 
the determination of the c-function is not straightforward 
and its computation far from fixed points is non-trivial 
even for very well known models, so that methods as form 
factor perturbation theory, truncated conformal space 
approach and conformal perturbation theory has been 
developed pQ. In d = 2 an expression of the c-function 
has been obtained in the framework of form factor per¬ 
turbation theory |7] for theories away from criticality and 
it has been applied to the sinh-Gordon model [5J. The 
sinh-Gordon model is a massive integrable scalar theory, 
with no phase transitions. In [8] one finds Ac = 1 for the 
sinh-Gordon theory. In a recent result [5], the analyti¬ 
cal continuation of the sinh-Gordon S'-matrix produces 
a roaming phenomenon exhibiting Ac = 1 and multi¬ 
ple plateaus of the c-function. The analytic continuation 
p —> i/3 of the sinh-Gordon model leads to the well-known 
sine-Gordon (SG) model with a periodic self interaction 
of the form cos (/ 3<t >). 

The SG model presents the unique feature to have 
a whole line of interacting fixed points with coupling 
(temperature) dependent critical exponents. It is in the 
same universality of the 2-dimensional Coulomb gas m 
and of the 2-dimensional AT EH, thus being one of 
the most relevant and studied 2-dimensional model, with 
applications ranging from the study of the Kosterlitz- 
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Thouless transition El to quantum integrability |12j and 
bosonization Il3l . In particular, for the SG model an 
ubiquitous issue is how to deal with the issue of the peri¬ 
odicity of the field El, which unveils and plays a crucial 
role for /3 7 ^ 0. Given the importance of the SG as a 
paradigmatic 2 -dimensional model, the determination of 
the c-function from the non-perturbative RG flow is a 
challenging goal, in particular to clarify the role played 
by the periodicity of the field for 0 ^ 0 . 

From the RG point of view, the determination of the 
behavior of the c-function is a challenging task requiring 
a general non-perturbative knowledge of the RG flow. 
Recently El. an expression for the Zamolodchikov’s c- 
function has been derived for 2 -dimensional models in 
the Functional RG (FRG) framework jlfiHTH] : resorting 
to an approximation well established and studied in the 
FRG, the Local Potential Approximation (LPA), an ap¬ 
proximated and concretely computable RG flow equation 
for the c-function was also written down |15j . By using 
this expression known results were recovered for scalar 
models on some special trajectories of the Ising and SG 
models. For the SG model, having a Lagrangian propor¬ 
tional to cos (/3</>), the determination and the integration 
of the c-function was carried out for 0 = 0 as a massive 
deformation of the Gaussian fixed point [T5] . Motivated 
by these results, both for the Ising and SG models and 
for general 2 -dimensional models, it would be highly de¬ 
sirable to have a complete description of the c-function 
on general RG trajectories. 

In the present paper we present the first numerical cal¬ 
culation of the c-function on the whole RG flow phase 
diagram of the SG model. The goal is to determine the 
behavior of the c-function, and the presence of known re¬ 
sults (namely, Ac = 1 ) helps to assess the validity of our 
approach along the different flows. We also complete the 
description initiated in El moving to more complex tra¬ 
jectories and showing that these cases are not a straight¬ 
forward generalization of the known results. We finally 
discuss the dependence of these results on the approxi¬ 
mation scheme used to compute FRG equations. 


II. FUNCTIONAL RENORMALIZATION 
GROUP METHOD 

In this section, we briefly summarize the FRG ap¬ 
proach [MU]- Starting from the usual concepts of Wil¬ 
son’s renormalization group (RG) it is possible to de¬ 
rive an exact flow equation for the effective action of any 
quantum field theory. This flow equation is commonly 
written in the form 


kd k T k [(f \ = ^Tr 


hdkRk 

_ r l 2) M + ’ 


(i) 


where is the effective action and F). ’[<p\ denotes the 
second functional derivative of the effective action. The 
trace Tr stands for an integration over all the degrees of 


freedom of the field <fi, while R k is a regulator function 
depending on the mode of the field and on the running 
scale k. When the running scale goes to zero k —> 0 
the scale-dependent effective action Ffc = o [<^] is the exact 
effective action of the considered quantum field theory. 

Usually equation (JT]) is treated in momentum space, 
thus the trace stands for a momentum integration and 
the regulator R k is a smooth function which freezes all 
the modes with momentum smaller than the scale k. 

The exact FRG equation (|T]) stands for functionals, 
thus it is handled by truncations. Truncated RG flows 
depend on the choice of the regulator function Rk, i.e. on 
the renormalization scheme. Regulator functions have 
already been discussed in the literature by introducing 
their dimensionless form 

Rk(p) =p 2 r{y), y=p 2 /k 2 , ( 2 ) 

where r(y) is dimensionless. Various types of regulator 
functions can be chosen (an archetype of regulator func¬ 
tions m has been shown to take the forms of the regula¬ 
tors used so far by setting its parameters). In this work 
we are going to consider the following regulators 

r pow {y) = (3a) 


r op t(y) = Q - 1^0(1 - y), (3b) 

where the first is known as the power-law type | 20 ] and 
the second one is the Litim (or optimized) |2 1 ) type regu¬ 
lator. Let us note that the so called mass cutoff regulator, 
which is used in [15], is identical to r pow (y) with 6 = 1 . 

One of the commonly used systematic approximations 
is the truncated derivative expansion, where the action 
is expanded in powers of the derivative of the field El: 



Vk(p) + Z k {p)-{d^p ) 2 H- 


In LPA higher derivative terms are neglected and the 
wave-function renormalization is set equal to constant, 
i.e. Zk = 1. In this case ([I]) reduces to the partial dif¬ 
ferential equation for the dimensionless blocked potential 
(14 = k~ 2 V k ) which has the following form in 2 dimen¬ 
sions 


(2 + kd k )V k (y) 


1 r°° y 2 ^f 

±_ / dy -- . 

47 r Jo (1 + r)y + V£{<p) 


(4) 


III. THE c-FUNCTION IN THE FRAMEWORK 
OF FUNCTIONAL RENORMALIZATION GROUP 


An expression for the c-function in FRG was recently 
developed in E- In this section we are going to give the 
guidelines of this derivation, reviewing the main results 
used in the next sections. 
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Let us start considering an effective action F < 7 ] for 
a single field <f> in curved space, with metric g^. We can 
study the behavior of this effective action under trans¬ 
formation of the field and the metric: 

(j) e^ T <(> (5) 

9 > e 2T g^ (6) 

where d$ is the conformal weight of the field = 
_d- 2 +j; £ or a sca j ar geld) while the background metric 
g^v has always conformal weight 2. From the require¬ 
ment that the effective action must be invariant under 
the Weyl transformation we obtain the following 

expression for a conformal field theory (CFT) in curved 
space na. 


a conformal invariant theory) and x is the fluctuation 
field. The notation T>Xd.b. stands for an integration over 
the fluctuation field x i n the curved space of the dila- 
ton background m- We can further manipulate latter 
expression moving cjjv on the l.h.s 

g — Sk [0,e 2 ' r 5] + (cE/v — Ck)Sp [e 2T S] _ 

T>Xd b e~ Suv ^ >+x ’ e ' 2Ts ^~^ Sk ^ x ’ e2T ^ ^ 

The Polyakov action in the dilaton background case as¬ 
sumes the form 

S p[ 9] = ~^ij tAt > ( 12 ) 


r[ <l>,g\ = ScFT[<t>,g\ + cSp[g]. (7) 

Scft[ 4>, g] is the curved space generalization of the stan¬ 
dard CFT action, which is recovered in the flat space case 
g ^„ = , c is the central charge of our theory and Sp [g\ 
is the Polyakov action term which is necessary to main¬ 
tain the Weyl invariance of the effective average action 
in curved space. 

To obtain FRG equations one has to add an infra-red 
(IR) cutoff term A S k [(j>, g] to the ultra-violet (UV) action 
of the theory. This is a mass term which depends both 
on the momentum of the excitations and on a cutoff scale 

k. 

A S k [<l>,g] = ^ J d 2 Xy/g<j)(x)R k {A)(j>(x), ( 8 ) 


where r is the dilaton field, A is the laplacian opera¬ 
tor and the integral is over an implicit spatial variable. 
Substituting latter expression into we obtain 




J Rxdb e~ SuV ^ +X,e2T ^~ ASk ^ X ’ e2 ^ S ] 


(13) 


In order to recover the usual flat metric integration we 
have to pursue a Weyl transformation ([5]) for the fields (j) 
and x 


e -S fc [e^ T 0,e 2T «5]- 


( c L/y- c fc) 

24 


f tAt _ 


J p^ e -S C /v[e^ T (0+x),e 2T 5]-AS fc [e^ T x ,e 2T 5] ) 


(14) 


where A is the spatial Laplacian operator. The effect is 
to freeze the excitation of momentum q <C k, but leaving 
the excitation at q > k almost untouched. The result of 
this modification of the UV action is to generate, after 
integrating over the field variable, a scale-dependent ef¬ 
fective action T k [<f>, g] which describes our theory at scale 
k. When the scale k is sent to zero the cutoff term in the 
UV action vanishes and the T k [(f>,g\ is the exact effective 
average action of the theory. 

The generalization of (J?]) in presence of the cutoff terms 
is 


and now the integration measure is in flat space. 

Finally deriving previous expression with respect to 
the logarithm of the FRG scale we obtain that the flow 
of the c-function dt.c k can be extracted from the flow of 
the cutoff action (| 8 ]) by taking the coefficient of the f r At 
term, 


kd k c k = 24TT(kd k AS k [e d,l,T x,e 2T 6}) 


f t At’ 


which after some manipulation becomes |15j 


(15) 


r fc [0, g\ = Sk [0, g\ + c k Sp [g\ H-, (9) 

where c k is now the scale-dependent c-function and the 
dots stands from some geometrical terms which do not 
depend on the field. We should now consider the case of a 
flat metric with a dilaton background g^ = e 2r Us¬ 
ing the standard path integral formalism for the effective 
action we can write 

e -r fc [0,e 2T <5] _ e -S k [^,e 2 ^8]-c k Sp[e 2T 5] _ 

T^Xd, b e ~ Suv ^~^ x,e2 ^^~ c vvSp[e 2 ' r 5]-ASfc[x,e 2 ' r 5] ^ ^ 

where Spv [<A> <?] is some UV action, cjjv is the value of 
the c-function in the UV (which can be equal to the cen¬ 
tral charge of some CFT if we are starting the flow from 


kd k c k = -\2-Kkd k G k {T] 


f T At 


(16) 


Eq. (16) shows that the c-function flow is proportional 
to the coefficient of the f tAt term in the expansion of 
the propagator flow kd k G k [T], also this flow has to be 
computed taking into account only the k dependence of 
the regulator function, i.e. 


kd k = kd k Rk 


d 

siv 


(17) 


This equation describes the exact flow of the c-function 
into the FRG framework. 

Since it is not in general possible to solve exactly equa¬ 
tion 0 and also equation (16) needs to be projected into 
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a simplified theory space to be computed numerically. In 
Ref. 15 j an explicit expression for the flow equation of 
the c-function in the LPA scheme has been derived with 
the mass cutoff 


kd k Ck = 


[kdkVk (Ao.fc )] 2 
(<Po,fc )] 5 


ya 


(18) 


with the dimensionless blocked potential V k (p) which is 
evaluated at its running minimum p = p 0ik (i.e. the 
solution of V{.{ip) = 0). We observe that an explicit ex¬ 
pression for the c-function beyond LPA is not available 
in literature. 

It should be noticed that, while Q is valid for any reg¬ 
ulator (cutoff) function, the expression for the c-function 
( fl8| ) has been obtained by using the mass cutoff, i.e. (3a) 
with 6=1. Other cutoff choices proved to be appar¬ 
ently very difficult to investigate. In the following, we 


will argue that while the expression (18) is sufficient to 


obtain a qualitative (and almost quantitative) picture of 
the c-function phase diagram the usage of other regulator 
functions is necessary to achieve full consistency. Where 
it is possible we will check the cutoff dependence of our 
numerical results. 


where /3 is scale-independent. 

In the mass cutoff case, i.e. the power law regulator 


(3a I with 6=1, one can derive [24] the flow equation for 
the Fourier amplitude of (21) from Eq. 


(2 + kd k )u k = 


1 


2TT/3 2 U k 


1 — a/1 — 


3±u 2 


( 22 ) 


(see Eq. (21) of [M] for vanishing mass). Similarly, using 
the optimized regulator (3b I gives 


(2 + kd k )u k = 


2n/3 2 u k 


1 


V 1 - < 


- 1 


(23) 


(19) 


B. The FRG equation for the SG model for 
scale-dependent frequency. 

A very simple, but still sensible, modification to ansatz 
is the inclusion of a scale dependent frequency, 


which, in order to explicitly preserve periodicity, should 
be rather considered as a running wave-function renor¬ 
malization. The ansatz then becomes 


IV. RG STUDY OF THE SINE-GORDON 
MODEL 


T k = d 2 x 


2 z k{dfj,<Px) 


Pfc (px 


(24) 


The SG scalar field theory is defined by the Euclidean 
action for d = 2 


r fc [<£] = J d 2 x 


7 ^{dij,<Px) 2 - ucos(/3p x ) 


(19) 


where f) and u are the dimensional couplings. Since we 
are interested in the FRG study of the SG model which is 
periodic in the field variable, the symmetry of the action 
under the transformation l22l 


p(x) -A p(x) + A 


( 20 ) 


is to be preserved by the blocking and the potential V k (p) 
must be periodic with period length A. It is actually 
obvious that the blocking, i.e. the transformation given 
by replacing the derivative with respect to the scale k by 
a finite difference in Q preserves the periodicity of the 
potential [221 [23] ■ 


A. The FRG equation for the SG model for 
scale-independent frequency. 

In LPA one should look for the solution of ([4]) among 
the periodic function which requires the use of a Fourier 
expansion. When considering a single Fourier mode, the 
scale-dependent blocked potential reads 


V k (p) = -u k cos (j3p), 


where the local potential contains a single Fourier mode 


V k {p) = -u k cos(p), 


(25) 


and the following notation has been introduced 

z k = 1//3 1 (26) 


via the rescaling of the field p —> <p/f3 k in (19), where z k 
plays the role of a field-independent wave-function renor¬ 
malization. Then Eq. ([I]) leads to the evolution equations 


kd k \4 = \ f p V k kd k R k 
kd k z k = T 0 V'" 2 J p V 2 kd k R k ( 


c) 2 -Dk „2 
dp 2 dp 2 F 




(27) 

(28) 


with V k = 1/ (z k p 2 +R k + V k ) and Vo = (27t)^ 1 f 0 n dip is 
the projection onto the field-independent subspace. The 
scale k covers the momentum interval from the UV cutoff 


A to zero. It is important to stress that Eqs. (27)-(28) are 
directly obtained using power-law cutoff functions. One 
may expect that these equations continue to be valid for 
a general cutoff provided that R k -A z k R k PTB]. This 
substitution has been tested for O(N) models, but its 
validity has been not yet discussed in the literature for 
the SG model. 

Inserting the ansatz (251 into Eqs. (27) and ( [28| the 
(21) RG flow equations for the coupling constants can be writ- 
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ten as 


kd k Uk = — 


1 

27r 


kOkRk 


Pk 


P 


Uk 




- 1 


(29) 


kd k z k = 


1 

27T 


I kd k R k i 

)p \ 


(ui P Hd p2 p k n4PZ+ui) 


4 (Pt-ulY 12 

u 2 k P k {d p 2 P k + p 2 d 2 2 P k )' 

2 ( P fc 2 -^) 5 / 2 


(30) 


with P k = + i?fc. In general, the momentum inte¬ 

grals have to be performed numerically, however in some 
cases analytical results are available. Indeed, by using 
the power-law regulator (3a I with 6=1, the momentum 
integrals can be performed [231 and the RG flow equa¬ 
tions read as 


(2 + kd k )u k = 
kd k z k = 


1 


27 TZ k U k 

1 


1 - 




24t r h 


(31) 


with the dimensionless coupling u = k 2 u. By using the 
replacements 


z k 


u k 


VPl 

Plu k , 


(32a) 

(32b) 


and keeping the frequency scale-independent ( dkZ k = 0 
i.e. d k /3 2 = 0) one recovers the corresponding LPA Eq. 


( 221 . 


V. c-FUNCTION OF THE SINE-GORDON 
MODEL FOR /3 = 0 


In this section we discuss the case /? = 0. We start 
by summarizing the results obtained for the c-function of 
the SG model in m The ansatz considered in m is 


mt 


VfcO) = -~B2 ( co s{Pk<p) - 1), 

Pk 


(33) 


where the frequency /3 k is assumed to be scale-dependent. 
If one directly substitutes (331 into the RG Eq. Q, 
then the l.h.s. of (4| generates non-periodic terms due 
to the scale-dependence of p k . Thus, the periodicity of 
the model is not preserved and one can use the Taylor 
expansion of the original periodic model 


Vk(<p) ~ - |y(wfe/3fc)</- 


(34) 


In this case, (33) is treated as a truncated Ising model 
and the RG equations for the coupling constants read as 


kdkrh k = 


kdkfil 


™l\Pk ~ 87r (! + **)] 


4tt( 1 + m£) 


1 (1 + 


l) 

4 fhl)Pt 


47 t (1 + m 2 ) 2 


(35) 

(36) 


The disadvantage of the scale-dependent frequency is 
that the periodicity of the model is violated changing 
the known phase structure of the SG model. However, 
the authors of m were interested in the massive defor¬ 
mation of the Gaussian fixed point which is at /3 = 0 and 
u = 0, so one has to take the limit (3 0 where the Tay¬ 

lor expansion represents a good approximation for the 
original SG model. Indeed, in the limit /3 —> 0, the RG 
Eqs. ([351) , (|36|) reduce to 


kdkfhl 

kdkPl 


™klPl -8^(4+ ™fe)] 


4tt(1 + rh 2 k ) 


— 2 m 2 


0 . 


fc (37) 
(38) 


Similar flow equations for the couplings ml and /3 k were 
given in m- The solution for the c-function based on 
(33) is in agreement with the known exact result, i.e. 


at the Gaussian UV fixed point cuv = 1 and in the IR 
limit cir = 0 , thus the exact result in case of the massive 
deformation of the Gaussian fixed point is Ac = 1 (Ac = 
cuv~cm). The numerical solution [15] gives Ac = 0.998 
in almost perfect agreement with the exact result. 

Although the numerical result obtained for the c- 
function in m is more than satisfactory, due to the Tay¬ 
lor expansion, the SG theory is considered as an Ising- 
type model. Thus, the RG study of the c-function start¬ 
ing from the Gaussian fixed point in the Taylor expanded 
SG model is essentially the same as that of the deforma¬ 
tion of the Ising Gaussian fixed point. So, it does not 
represent an independent check of (|18[ ). Indeed, insert¬ 
ing (371 into (18) using the ansatz (1331) one finds 


kd k c k = 


4 m| 




(39) 


which is identical to Eq. (5.3) of [15] (with a = 1) ob¬ 
tained for the massive deformation of the Gaussian fixed 
point in the Ising model and it can be also derived from 
Eq. (5.19) of |TS] in the limit of /3 2 —> 0. 

Therefore, it is a relevant question whether one can re¬ 
produce the numerical results obtained for the c-function 
(with the same accuracy) if the SG model is treated 
with scale-independent frequency (21), or beyond LPA, 


by the rescaling of the field (25). Also ref. [15] treats 
only massive deformations of non interacting UV fixed 
points, then on such trajectories only the mass coupling 
is running. Nevertheless the c-theorem should hold on 
all trajectories, even when more couplings are present. 
Our aim is to demonstrate that the derivation of m is 
valid even in these more general cases, but, due to trun¬ 
cation approach, the approximated FRG phase diagram 
does not fulfill the requirements of the c-theorem exactly 
and, therefore, only approximated results are possible. 


VI. c-FUNCTION OF THE SINE-GORDON 
MODEL ON THE WHOLE FLOW DIAGRAM 


In this section we study the c-function of the SG model 
on the whole phase diagram, studying both the scale in- 
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dependent wave-function renormalization and the treat¬ 
ment with the running frequency. 


A. Scale-independent frequency case 


The definition for the SG model used in this work, i.e. 


(19), differs from (33) because the frequency parameter 


is assumed to be scale-independent in LPA. The running 
of p can only be achieved beyond LPA by incorporating 
a wave-function renormalization and using a rescaling of 
the field variable which gives Zk = l//3f. 

Let us first discuss the results of LPA. Equations (|22|) 


and (23) have the same qualitative solution. In Fig.lllwe 


show the phase structure obtained by solving (23). 



k/A 



(f/ 7T 


FIG. 1: The figure shows the phase structure of the SG 
model obtained by the FRG equation using Litim’s regulator 
in the scale-independent frequency case. The two phases are 
separated by p 2 = 87r. The dashed line shows the line of IR 
fixed points of the broken phase. 


The RG trajectories are straight lines because in LPA 
the frequency parameter of (19) is scale independent. 
Above (below) the critical frequency /3 2 = 87r, the line 
of IR fixed points is at uir = 0 (uir ^ 0). For /3 2 < 87r 
the IR value for the Fourier amplitude depends on the 
particular value of /3 2 thus, one finds different IR effec¬ 
tive theories, i.e. the corresponding CFT depends on the 
frequency too. 

The scaling for the c-function is the one expected from 
the c-theorem. It is a decreasing function of the scale 
k which is constant in the UV and IR limits, see Fig.jT] 
Due to the approximation of scale independent frequency 
P, here, the IR value of the c-function depends on the 
particular initial condition for /3 2 . Then when we start at 
the Gaussian fixed points line (c = 1), in the symmetric 
phase, the flow evolves towards an IR fixed point, but at 
this approximation level, we have different IR fixed points 
which are all at different u values and consequently the 


FIG. 2: Running of the c-function obtained in the scale- 
independent freq uenc y case by solving (22 I (dotted lines - 
mass cutoff) and (231 (dot-dashed lines - Litim cutoff) com¬ 


bined with (|18| for the SG model is plotted for various val¬ 
ues of the frequency ft 2 . From top to bottom it is P 2 /n = 
4,2,0.005. Due to the poor convergence properties of ( |22| , 
where the mass cutoff was used, the RG flow stops at some 
finite momentum scale and the deep IR value of the c-function 
cannot be reached (dotted lines). The use of the Litim cutoff 
( |23[ ) (dot-dashed lines) can produce us the IR constant for 
the c-function. However, for very small value of /3 2 the low- 
frequency approximation is the best choice, i.e. one has to 
solve ( |41| (green line): results from different cutoff functions 
are indistinguishable. The solid lines represent the results ob¬ 
tained with power law cutoff (6 = 2), while the dashed lines 
are the results with the exponential cutoff. The inset shows 
the results for an enlarged theory space where higher harmon¬ 
ics are included in (|25|) (dot-dashed lines - Litim cutoff). 


Ac values differ from the exact one. The exact result 
Ac ~ 1 is obtained only in the p —>■ 0 limit. 

We notice that Eq. (22), where the mass cutoff was 
used, has very poor convergence properties and the flow, 
obtained from them, stops at some finite scale, thus 
the deep IR values of the c-function cannot be reached 
(dashed lines in Fig. [2]). 

The use of the Litim cutoff RG Eq. (231 improves the 
convergence of the RG flow but the IR results for the 
c-function are very far from the expected Ac = 1, which 
can be recovered only in the vanishing frequency limit. 
Also the inclusion of higher harmonics in (25) ( inset in 
Fig{2]) does not improve this result. 

It should be noted that Eq. (18) is strictly valid only 
in the mass cutoff case, however in Fig. [2] we used Eq. 
(18) even in the optimized cutoff case. This inconsistence 
cannot be regarded as the cause for the unsatisfactory 
results obtained in the large P cases, indeed we expect 
very small dependence of the flow trajectories upon the 
cutoff choice. 

This small dependence on the regulator is evident from 
the comparison of the mass and Litim regulator results 
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of trajectories for the c-function in Fig. [2j which are 
very similar, at least in the region where no convergence 
problems are found. This similarity justifies the use of the 
mass cutoff result ( fl8] ) with RG flow Eqs. (23) obtained 
by the optimized (Litim) regulator. 

We also computed the c-function flow for other cutoff 
functions, namely the power-law 6 = 2 (solid lines in Fig. 
[ 2 ]) and the exponential one (dashed lines). Apart from 
the mass cutoff, all the others converge to the IR fixed 
point. The conclusion is that there is not a pronounced 
dependence of the findings on the cutoff schemes and that 
the constant frequency case in not sufficient to recover 
the correct behavior for the c-function. 

We observe that the lack of convergence observed in 
mass cutoff case is not present in the small frequency 
limit analyzed in m- Indeed, expanding flow equations 
( 22 ) and (23) we get 


kd k Uk ~ - 2 iik + 


u k /3 2 

4-7T 


2 Uk ; 


(40) 


which is valid for vanishing frequency and it is indepen¬ 
dent of the particular choice of the regulator function, 
i.e. it is the same for the mass and Litim cutoffs. Sub¬ 


stituting (40) into (18) using (21) the following equation 


is obtained for the c-function of the SG model: 


kOkCk 


(kdkUk/3 2 


(-2 u k /3 2 




(1 + u k /3 2 ) {l + Uk/3 2 ) 


m 


213 


(41) 


where the identification rh 2 = UkP 2 is used. The scale 
dependence of the c-function in that case is identical to 
the massive deformation of the Gaussian fixed point and 
the corresponding RG trajectory is indicated by the green 
line in Fig. [3] 

It is important to note that for finite frequencies /3 2 ^ 0 


the Taylor expanded potential (34) cannot be used to 


determine the c-function since it violates the periodicity 
of the model. In this case only Eqs. (22) or (23) can 
produce reliable results. 

In order to improve the LPA result for the c-function 
of the SG model without violating the periodicity of 
the model one has to incorporate a scale-dependent fre¬ 
quency, i.e. a wave-function renormalization (we refer to 
this approximation as 2 +LPA), as it is discussed in the 
next subsection. 


B. The scale-dependent wave function 
renormalization 

The inclusion of the running wave-function renormal¬ 
ization changes the whole picture of the SG phase dia¬ 
gram, with all the « / 0 fixed points collapsing into a 
single (/3fc = 0 , u = 1 ) fixed point, as it is expected from 
the exact CFT solution. The phase diagram obtained at 
this approximation level is sketched in FigjdJ where we 
evidence three different regions. The Ac is strictly well 
defined only in region /, where we start from a Gaussian 



FIG. 3: The flow diagram of the SG model in the scale- 
dependent frequency approximation. The phase space is di¬ 
vided into three regions. In region I we have a line of UV 
repulsive Gaussian fixed points (u = 0,/3 2 < 87r). Every tra¬ 
jectory starting in the vicinity of this line ends in an IR at¬ 
tractive fixed point (purple full circle, u = 1, /3 2 =0). The Ac 
observed along the trajectories of this region should be equal 
to 1. Region II contains a line of IR attractive Gaussian fixed 
points (u = 0, /3 2 > 87t,) which are the end points of trajec¬ 
tories starting at /3 2 « 00 below the thick green line, i.e. the 
separatrix. Region III contains those trajectories starting at 
P 2 ~ 00 which end in the IR attractive fixed point (purple 
full circle). 


fixed point Cuv = 1 and we end up on a massive IR fixed 
point cjr = 0. The massive IR fixed point related to the 
degeneracy of the blocked action is an important feature 
of the exact RG flow [26H29] and it was considered in SG 
type models [551 13 H 1131 ] ■ 

In region II the trajectories end in the Gaussian fixed 
points c = 1 but they are coming from infinity where 
actually no fixed point is present. This is due to the 
fact that we are not considering in our ansatz (24) any 
operator which can generate a fixed point at c > 1 and 
then the trajectories ending at c = 1 are forced to start 
at infinity. Thus, Ac is not defined in this region. 

Region III contains those trajectories which start at 
f3 = 00 but end in the IR massive fixed point at c = 0. 
Even in this case the Ac is not well defined. 

In the following we are going to discuss in details the 
results of region / where all the trajectories should give 
Ac = 1. We shall ignore region II where the Ac is not 
defined, briefly discussing region III. 

The presence of wave-function renormalization is nec¬ 
essary to obtain the qualitative correct flow diagram for 
the SG. Note that Eq. (18) has been derived only in the 
case of scale-independent kinetic term and the deriva¬ 
tion of an equivalent expression in the case of running 
wave function renormalization appears far more demand¬ 
ing than the calculation sketched in this paper. How¬ 
ever, it is still possible to get a sensible result using the 



























mapping between the running wave-function renormal¬ 
ization and the running frequency [3 k cases (as shown in 
Eqs. (32a) and (32b)), finally obtaining Eq. ( |44| In other 
words the equation (181 is valid only at LPA level, but it 


is still possible to apply it to the z+lA’A sch eme, since , 
thanks to the mapping described in Eqs. (32a) and (32b) 


the z+LPA ansatz can be mapped into an LPA one. 
We will then use directly the ansatz, 

V k = u k cos ((3 k cj)) 


(42) 


with no wave function renormalization present in the ki¬ 
netic term. This ansatz is equivalent to ansatz (241 if we 


rescale the field and use the relations (32a) and (32b), 


with the running frequency playing the role of a wave- 
function renormalization. 


Ansatz (42) is not suited to study the SG model when 


full periodicity has to be preserved, indeed when we sub¬ 
stitute it into equation (27) symmetry breaking terms 


appear. The same happens when we substitute it into 
Eq.(18). However in the latter case symmetry breaking 


terms are not dangerous, since we have to evaluate the 
expression at the potential minimum where all the sym¬ 
metry breaking terms vanish. 

Proceeding in this way we obtain 


kd k c k = 


(PlkdkUk + 2u k (3 k kd k (3 k ) 2 


(1 


' « fc / 3 fe ) 3 


(43) 


where no inconsistency is present. 


We still cannot use expression (43), since we cannot 


write a flow for (3 k due to the non-periodic terms. To 


avoid these difficulties we rewrite expression (43) using 
the inverse transformation of (|32a) and (32bI, 


kd k c k = 


( kd k u k ) 2 

(1 + Mfe) 3 ' 


(44) 


The last expression is fully coherent and represents the 
flow of the c-function in presence of a running wave- 
function renormalization into the SG model; it is worth 


noting that the use of transformations (32a) and (32b) 


gave us the possibility to derive the expression (43) from 


Eq (18), which was derived in m in the case of no-wave- 
function renormalization. 

In the limit fi‘ k=1 \ —> 0, the IR result of the c-function 
(see the inset of Fig. [4]) tends to zero. This implies that in 
the limit /3^_ A —► 0 the difference Ac —> 1. The numerical 
result found in this case reaches the accuracy 1 > Ac > 
0.99 of the scale-independent frequency solution (33) but 


now the periodicity of the SG model is fully preserved 
(which was not the case in )J5]). It should be also noted 
that the accurate results of F ig|4| could not be obtained in 
the mass cutoff framework ( |3a[ ) with 6 = 1, which does 
not allow the flow to converge, but our findings were 
obtained with the smoother 6 = 2 cutoff. 

Fig [4] reports the running of the c-function for various 
values of the initial condition (3\. The final Ac value 
depends on the trajectory even if it should not be at exact 



k/A 

FIG. 4: Running of the c-function obtained in the case of 
scale-dependent wave-function renormalization for the single¬ 
frequency SG model, as expected the case of small frequency 
(l3 k =A < 0.1) was already very well described by the scale- 
independent frequency case Fig.[T| The inset shows the results 
obtained for cir as a function of p k =A, these results show 
lower accuracy in the large frequency limit, while they become 
practically exact in the limit /3 k =A —> 0, accordingly with Di- 


level. This discrepancy shows that the flow obtained by 
approximated FRG procedure cannot satisfy the exact 
CFT requirements for the c-function. 

The discrepancy between the exact Ac = 1 value and 
the actual results obtained by the FRG approach can be 
used to quantify the error committed by the truncation 
ansatz in the description of the exact RG trajectories. 

We observe that the results of Fig. [4] main and inset 
are obtained by using power-law regulator with 6 = 2. 
The same computation appears to be considerably more 
difficult using general cutoff functions, including the ex¬ 
ponential one. 

Let us note that for vanishing frequency the RG flow 
equations become regulator-independent and that the c- 
function value tends to the exact result Ac = 1. This 
justifies the accuracy obtained in m even though the 
mass cutoff was used and the periodicity violated. 

Finally we go on showing the results in region III. 
As discussed in the description of Fig.[3j trajectories in 
region III of the SG flow diagram should not have a well 
defined value for the c-function, due to the fact that those 
trajectories start at (3 k= A = 00 where no real fixed point 
is present. 

However the numerical results obtained for those tra¬ 
jectories Fig. [5] are not so far from Ac = 1, due to the 
fact that they get most of the contribution in the region 
where they approach the ’’master trajectory” separatrix 
of region I i.e. the blue thick line in Figj3j which we 
know to have a value Ac ss 1, while the portion of the 
trajectories close to region II get almost zero contribute. 
The results of region III are also in agreement with the 
findings of region II (not shown) where Ac ~ 0 for all 
the initial conditions. 
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k/A 

FIG. 5: The flow of the c-function in Region III of the SG 
flow diagram Fig. [3] the result is approximately Ac = 1 due 
to the fact that in region III most of the contribution to the 
c-function comes from the part of the trajectories very close 
to ’’master trajectory” separatrix of region I (the blue thick 
line in Fig.[3|. 


VII. CONCLUSIONS 


In this paper we provided an estimation of the c- 
function over the RG trajectories of the sine-Gordon (SG) 
model in the whole parameter space. Using this result we 
showed that the numerical functional RG study of the 
SG model with scale-dependent frequency recovers for 
/ 3 2 < 8n (region I of Fig. til) the exact result Ac = 1 with 
a good quantitative agreement while preserving the pe¬ 
riodicity, which is the peculiar symmetry of the model. 
We also pointed out the dependence of this c-function 
calculation on the approximation level considered. For 
p = 0 one retrieves directly Ac = 1, also in the scale- 
independent frequency case, while for P ^ 0, again using 
scale independent frequency, we recover this result in the 
p 2 —>• 0 limit, while increasing P 2 up to 8tt in region I as a 
result of the used approximation the agreement becomes 
worst, remaining anyway reasonably good, as shown in 

Fig-0 

Retrieving Ac = 1 is the SG counterpart of the com¬ 
putation of Ac for the sinh-Gordon model [8j [9]. This 
result can be understood by noticing that the analytical 
continuation /? -+ ip [32] may be expected not to alter 
the Ac defined in the Zamolodcliikov theorem, and that 
functional RG even in its crudest approximation does not 
spoil such correspondence for Ac, provided that the pe¬ 
riodicity of the SG field is correctly taken into account. 

We developed a fully coherent expression for the c- 
function in the case of running frequency, which gives 
better results in the whole region I (defined in Fig. [3]). 
These results are compatible with the exact scenario up 
to an accuracy of 80% in the whole region I. Such accu¬ 
racy grows to 99% in the small beta region in agreement 
with [T5], as we discussed in Fig.[4j 


We also noticed that for P 2 > 0 the use of the mass 


cutoff, as necessary to be consistent with expression (181, 


is not possible due to bad convergence properties, and 
the use of different b values or of different cutoff types is 
needed. 

It should be noted that while the numerical results are 
quite accurate the exact property that all trajectories 
of region I should have the same value for Ac is not 
preserved by truncations schemes (24) (results in Fig. [l]) 
nor (19) (results in Fig.[4|. Actually the P 2 —> 0 limit 
always gives the correct result even when treated with 
the most rough truncation, this result being independent 
from the cutoff function. At variance one needs to go to 
the running frequency case to obtain reliable results for 
p 2 » 0. 

Even when full periodicity in the field is maintained, 
the z+LPA truncation scheme is not sufficient to recover 
exact results for the c-function. Indeed the quantity com¬ 
puted using expression (18) satisfies two requirements of 


Zamalodchikov’s c-theorem, 

1 . dtCt > 0 along the flow lines, 

2 . dtCt = 0 at the fixed points, 

but it fails to reproduce the exact central charge of SG 
theory. 

The final result of our calculation also depends on the 
chosen cutoff function and, as already mentioned, it was 
not possible to use the same cutoff scheme for both the 
couplings flow equations and the c-function flow (181. We 


do not expect these issues to be responsible for the error 
in the fixed point value of the c-function. Modifications 
of the cutoff scheme in LPA calculations have small influ¬ 
ence on the results (around 5%) and we may expect this 
property to be maintained at z+LPA level, where calcu¬ 
lations with different cutoff functions were not possible. 

The main source of deviation from the exact result 
Ac = 1 is then probably due to z+LPA truncation in 
itself. We are not able to identify whether this devia¬ 
tion is only due to the approximation in the c-function 
flow or rather to the description of the fixed point given 
in z+LPA, which does not reproduce the exact central 
charge. 

Certainly the c-function flow at z+LPA level not 
merely violates the exact fixed point value of the c- 
function, but it is also not able to produce trajectory 
independent results, as shown in Fig. [4] This scenario 
is not consistent with an unique central charge value at 
the SG fixed point and it is then impossible to dig out 
any information about this quantity from this approach. 
In this perspective it would be interesting to have an 
independent method to calculate the fixed point central 
charge at a given truncation level. 

Obviously the reproduction of the Zamalodchikov’s c- 
theorem should be better satisfied increasing the trunca¬ 
tion level considered. However it has been shown that, at 
LPA level, the addition of further harmonics in the po¬ 
tential does not improve the results presented, while the 
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introduction of running frequency is crucial to achieve 
consistency of the phase diagram. 

This situation is peculiar of the LPA truncation level. 
Beyond 2 +LPA we expect the most relevant corrections 
from higher harmonics in the potential and only small 
variations are expected from the introduction of higher 
fields derivatives in (24) 

Finally we remark that the trajectories of the other two 
regions do not have a definite Ac value, however while 
region II gives results Ac « 0, region III has the Ac 
values close to the ones obtained in region I Fig.[5j due 
to the fact that all the trajectories in this region merge 
with the ’’master trajectory” separatrix of region I in the 
k —> 0 limit. 

We conclude by observing that in our opinion a rel¬ 
evant future extension of this work could be the study 
of the c-function for the Ising model and for minimal 
conformal models in general. Our work also points out 
to the possible investigation with RG techniques of the 
analytical continuation relating the sine- and the sinh- 
Gordon models. From this respect we think it could be 


worthwhile to systematically study models interpolating 
between these two celebrated cases, both to highlight the 
roaming phenomenon for integrable interpolations and to 
put forward critical properties of non-integrable interpo¬ 
lations. 
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